REPRESENTATIONS OF NON-NEGATIVE POLYNOMIALS 

VIA CRITICAL IDEALS 



DANG TUAN HIEP 

Abstract. This paper studies the representations of a non-negative 
polynomial / on a non-compact semi-algebraic set K modulo its critical 
ideal. Under the assumptions that the semi- algebraic set K is regular 
and / satisfies the boundary Hessian conditions (BHC) at each zero of 
/ in K, we show that / can be represented as a sum of squares (SOS) of 
real polynomials modulo its critical ideal if / > on if . In particular, 
we focus on the polynomial ring M\x\. 



1. INTRODUCTION 

We know that a polynomial in one variable f{x) G M[x] satisfies f{x) > 0, 
for all X G M, then /(x) = Xll^i fl'i^(^)' where gi{x) G M[x], i.e., / is a sum 
of squares in M.[x] (SOS for short). However, in multi-variate cases, this is 
not true. A counterexample was given by Motzkin in 1967. If f{x,y) = 
1 -|- x'^y^ + x^y* — Sx'^y^, then f{x,y) > 0, for all x,y E M.. But / is not 
a SOS in M[x,?/]. To remedy that, we will consider the polynomials that 
are positive on K, where i^' is a semi-algebraic set in M". For example, 
Schmiidgen's theorem [Schmj states that for a compact semi-algebraic set, 
every strictly positive polynomial belongs to the corresponding finitely gen- 
erated preordering. Afterward, Putinar |Pu] simplified this representation 
under an additional assumption by using the quadratic module instead of 
the preordering. However, these results of Schmiidgen and Putinar have two 
restrictions. Firstly, the polynomials are positive, not merely non-negative. 
Secondly, K must be a compact semi- algebraic set. Hence we seek to iden- 
tify the representations of the non-negative polynomials on the non-compact 
semi-algebraic sets. 

In |NDSj . the authors presented a representation of the non-negative 
polynomials on the whole space modulo their gradient ideals. Afterward, in 
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[DNPj . the authors proved a similar representation on the arbitrary semi- 
algebraic sets. These results were achieved under the condition of the cor- 
responding ideals must be radical. However, it is not simple to check this 
condition. In order to overcome such limitation, in jM], Marshall considered 
another condition - the boundary Hessian condition (BHC). He proved that 
the result in |NDSj still held true if the radical condition is replaced by the 
BHC condition. In [Hi], the author presented an extension of theorem 2.1 
in [M] in the same way that the result in |DNPj was the extension of the 
corresponding result in |NDS] . 

However, in [Hi] and jDNPj the authors considered a larger polynomial 
ring ]R[x, A], i.e., they added Lagrange multipliers to the representations. 
This paper will help us overcome this. We will present the representations 
of the non-negative polynomials via their critical ideals. In particular, we 
focus on the polynomial ring M.[x]. 

2. Preliminaries 

In this section, we present some notions and results from algebraic geom- 
etry and real algebra needed for our discussions. The readers may consult 
[BUR] . [ULU] , and [PD] for more details. 

Throughout this paper, denote by M.[x] the ring of polynomials in x = 
{xi, . . . ,Xm) with real coefficients. Given an ideal / C M.[x], define its com- 
plex variety to be the set 

V{I) = {x G C" I p{x) = 0, Vp G /}, 

and its real variety to be 

v^{i) = v{i) n M™. 

A nonempty variety V = V{I) C is irreducible if there do not exist 
two proper subvarieties Vi,V2 C V such that V = Vi U V2. The readers 
should note that in this paper, "irreducible" means that the set of complex 
zeros cannot be written as a proper union of subvarieties defined by real 
polynomials. 

Given any ideal / of M[a;], its radical ideal vT is defined to be the fol- 
lowing ideal: 

a/7 = {g G R[x] \ q'' e I for some / G N}. 
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Clearly, I C y/l; I is a radical ideal if a/7 = I. As usual, for a variety 
V C C"^, I{V) denotes the ideal in C[x] of polynomials vanishing on V. We 
will write I*{V) for the ideal I{V) n R[x]. 

We need versions of the NuUstellensatz for varieties defined by polyno- 
mials in M[a;]. The following two theorems are normally stated for ideals in 
C[a;]; however, keeping in mind that V{I) lies in C™", they hold as stated 
for ideals in M.[x]. 

Theorem 2.1 f fCLO] ). // / is an ideal m R[x] such that V{I) = 0, then 

1 e /. 

Theorem 2.2 f fCLQ] ). /// zs an ideal m R[x], then I^{V{I)) = Vl. 

Let Qi, . . . ,gs G M.[x]. We define the preordering generated hj gi, . . . , Qs 
as follows: 

^=1 E ^eg1'---9T], 

where e = (ei, . . . , e^) G {0, 1}* and CTg are sums of squares of polynomials 
in R[x]. 

We also define the semi-algebraic set generated hj gi, . . . , gs as follows: 
= {x G M" I giix) > 0, i = 1, . . . , s}. 

Definition 2.1 (see [NW] . Definition 12.1). For each x G M", let be 
the set of indices j for which gj vanishes at x. The semi-algebraic set K is 
called regular, if for each x E K, the vectors Vgj{x),j G Jx, are linearly 
independent. 

Throughout this paper, we always assume that the semi-algebraic set K 
is regular. 

3. The critical variety 
Definition 3.1. The critical variety of f on K is defined as follows: 
C{f,K) := {x G M" I there exist real numbers Aj such that 

s 

Vf{x)-J2^^^9^{x)=0, 

i=l 

\9iix) = 0,2 = 1, . . . ,s}. 

Remark 3.1. 
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(i) In the global case, i.e., when the semi-algebraic set K is the whole 
space M", we have 

C(/,K) = {xGR'M V/(x)=0}, 

which is the real gradient variety of / (see |NDSj ). 

(ii) Consider the projection tt: M" x — )■ M", (x. A) i— t- x, where vari- 
ables A = (Ai,...,As) are Lagrange multipliers. Then C{f,K) = 
vr(VKKT), here 

s 

^KKT:={(x,A)eM"xM^ I V/(x)-^AiV(?.(x) = 0, 

1=1 

XiQiix) = 0,2 = 1, . . . ,S}, 

is the real KKT variety of f on K (see |DNPj ). 
In this section, we will study the properties of the critical variety C(/, K). 

Proposition 3.1. The following statements hold true 

(i) C(/, K) = C{f + a, K), for all a G M. 

(ii) /// attains its infimum at x* G K, then x* G C{f,K). 

Proof. 

(i) We see clearly that V/ = V(/ + a), for all a G M. Then, by definition 
of the critical variety, we have C{f, K) = C{f + a, K), for all a G M. 

(ii) By Karush-Kuhn- Tucker theorem (see e.g. |NWj ) . if / attains its 
infimum at x* G K, then there exist Ag, A]^, . . . , A* at least one of 
which is different from zero, such that 

s 

KVf{x*)-Y,X*Vg,{x*) = 0, 

i=l 

Since K is regular, then we can choose Aq = 1. 
Thus X* G C{f,K). 

□ 

We will use the following notations in the remainder of the paper. 
Definition 3.2. For each subset J of {1, ... , s}, we consider the polynomial 

«.fx)-l n.^^^^-^"^) ' 
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If J = {ji ■ . . ,jk}, we will denote by hj G M[x] the following polynomial 

hj{x) := det{Aj{x)A^{x)), 

where 



Aj(x) :-- 



( dj_ dL ... 9f \ 

dxi dx2 dxn 

ddh ^9h . . . 

dxi dx2 dxn 



\ dxi dx2 dxn / 



zs a (fc + 1) X n-matrix. Observe that hj{x) = if and only if the vectors 
VfyVgjjj G J are linearly dependent. 

Proposition 3.2. The critical variety C{f,K) is an algebraic set. More 
precisely we have 

C{f, K) = {xeR^\ gj{x)hAx) = 0, VJ C {1, . . . , s}}, 

where we use the notation := {1, . . . , s}\J. 

Proof. The proof is similar as that of Proposition 3.1 in |HP2] and therefore 
is omitted here. □ 

4. Boundary Hessian Conditions, gradient ideals and KKT 

IDEALS 

We say / satisfies the BHC (boundary Hessian conditions) at the point 
X* in K if there are some k E {1, . . . , n}, and vi, Vk E N with 1 < vi < 
... < Vk < s such that g^^^, . . . , g^^ are parts of a system of local parameters 
at X*, and the standard sufficient conditions for a local minimum of at x* 
hold, where L is the subset of defined hj g^^^^x) > 0, . . . , g^i, (x) > 0. This 
means that if ti, . . . , t„ are local parameters at x* chosen so that tj = g^. for 
i < k, then in the completion . . . , of M[x] at x*, / decomposes as 

/ = /o + /i + /2 H — ■ (where fi is homogeneous of degree i in the variables 
ti, . . . ,tn with coefficients in M), /i = aiti + ■ ■ ■ + aktk with > 0,i = 
1, . . . ,k, and the (n— A;)-dimensional quadratic form /2(0, . . . , 0, t^+i, . . . , t^) 
is positive definite. 

Theorem 4.1 (Marshall [M])- If f satisfies the BHC at each zero of f in 
K, then f G P + 

Example 4.1. Let f,gi G ]R[x,?/, z] be given by 

/(a;, y, z) = x; gi{x, y,z) = x - y'^ - . 
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Then 



K = {{x,y, z) E \ z - - z'^ > 0}. 



Clearly, / > on i^T, and the unique zero of / in K occurs at (0,0,0). 
Furthermore, / satisfies the BHC at (0, 0, 0). Indeed, let ti = gi = x — y"^ — 
z^i t2 = y and ^3 = z. These form a system of local parameters at (0, 0, 0). 
Then f = x = {x -y^ - z^) +y^ + z^, so f = /1 + /2, where fi{ti,t2,h) = h, 
and f2(ti,t2,ts) = ^2 + ''^i- Also, the coefficient of ti in fi is positive (it is 
1), and ^2,^3 do not appear in fi. The quadratic form 72(0,^25 ^3) = ^2 + ^3 
is positive definite (when viewed as a quadratic form in the two variables 
^2,^3)- So, according to the definition, / satisfies the BHC at (0,0,0). Here 
/ has a representation as follows: 

/ = ctq + 0-15(1 + hf, 

where ctq = y"^ + z'^, ai = 1, h = 0. 

Now we define the gradient ideal of / as follows: 



Under the assumption that Jgrad is radical, we have the following result. 

Theorem 4.2 (Nie-Demmel-Sturmfels |NDSj ). Suppose that 
(i) f >0 on M", 
(a) Jgrad is radical. 
Then f is a sum of squares modulo Jgrad- 

If we replace the radical condition of Jgrad by an another condition that 
/ satisfies the BHC at each zero of /, then we will have the following result. 

Theorem 4.3 (Marshall [M])- Suppose that 
(i) f >0 on R", 

(a) f satisfies the BHC at each zero of f . 
Then f is a sum of squares modulo Jgrad- 

Similar to generalization of the gradient ideal, we define the KKT ideal 
of / as follows: 




Jkkt — (J^i, • • • , 



where 
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Two following results are generalizations of theorem 14.21 and theorem 14.31 
in the same way. 

Theorem 4.4 (Demmel-Nie-Powers |DNPj ). Suppose that 
(i) f >0 on K, 
(a) /kkt is radical. 
Then f E P + /kkt- 

Theorem 4.5 (Hiep [Hi]). Suppose that 
(i) f >0 on K, 

(a) f satisfies the BHC at each zero of f in K. 
Then feP + /kkt. 

Remark 4.1. The radical condition and the BHC condition are different. 
This means that there exist polynomials which satisfy the radical condition, 
but do not satisfy the BHC condition and conversely. The following example 
will demonstrate this difference. 

Example 4.2 (Marshall [M]). 

1. Let n = 1 and s = (so that K = M.). Then the polynomial in one 
variable f{x) = 6x^ + 8a;'^ + 3a;^ satisfies the BHC condition, but 

it does not satisfy the radical condition. Indeed, — = 12x{x + 1)^, 

f(x) > on M, / has a zero at x = 0, and ^—^(0) = 12 > 0. However, 
the gradient ideal / = (12a;(a; + 1)^) which also is the KKT ideal, is 
not radical, because g{x) = x{x + 1) G vT, but g ^ I. 

2. Let n = 2 and s = (so that K = R^). Then the polynomial in 
two variables f{x,y) = x^ does not satisfy the BHC condition, but 
it satisfies the radical condition. Indeed, the Hessian matrix of / is 
not positive definite at any zero of / in K. However, the gradient 
ideal / = (2x) which also is the KKT ideal, is radical. 

Remark 4.2. If we leave both the radical condition and the BHC condi- 
tion, then we will have the corresponding representations of strictly positive 
polynomials. 

Theorem 4.6 (Nie-Demmel-Sturmfels [NDS]). // / > on W, then f is a 
sum of squares modulo /grad- 

Theorem 4.7 (Demmel-Nie-Powers |DNP] ). If f > on K , then f G 
P + /kkt ■ 
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Remark 4.3. In the proof of theorem 14.41 theorem 14.51 and theorem 14.71 
we must work in a larger polynomial ring R[x, A], i.e., we must add the 
Lagrange multipliers to our representations. 

5. Sums of squares modulo critical ideals 

In this section, we present our main results. These are similar to theorem 
14.41 and theorem 14. 5[ but without modulo /kkt- It is replaced by modulo 
another ideal - the critical ideal of / on K. In its proof, we work particularly 
in the polynomial ring M[a;]. 

Let us start with some notations. The ideal 

lU.K) := ((7>jc,VJC{l,...,s}) 

generated by gjhjc is called the critical ideal of f on K. By Proposition 
13. 2^ we have 

CU.K) = V^{IU.K)). 

Theorem 5.1. Suppose that 
(i) f>0 on K, 

(a) f satisfies the BHC at each zero of f in K. 
Then f eP + I{f, K). 

To prove the theorem 15. we need the following lemma. 

Lemma 5.1. Let W he an irreducible component ofV{I{f, K)). //W^flM" 7^ 
0, then f is constant on W . 

Proof. This follows from the proof of lemma 3.6 in [HP3j . □ 

Proof of theorem \5. 1\ We decompose V{I{f, K)) into its irreducible compo- 
nents and let Wq be the union of all the components whose intersection with 
K is empty. We note that this includes all components W with ly flR" = 0. 
Thus, by lemma WA\ f is constant on each of the remaining components. We 
group together all components for which / takes the same value. Then we 
have pairwise-disjoint subsets Wi, . . . , Wr of W such that for each i, / takes 
a constant value on Wi, with the Oj being distinct. Further, since each 
Wi contains a real point and / is non- negative on C(/, K) fl the value of 
/ on each Wi is real and non-negative. We assume ai > ■ ■ ■ > > 0. We 
fix a primary decomposition of /(/, i^), for each i G {0, 1..., r}, let Jj be the 
intersection of those primary components corresponding to the irreducible 
components occurring in Wi. Thus, V{.Ji) = Wj, Vi = 0, 1, . . . , r. 
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Since Wi H Wj = 0, we have Jj + Jj = M.[x] by theorem 12.11 There- 
fore the Chinese remainder theorem (see, e.g., jE]) imphes that there is an 
isomorphism 

ip : R[x]/I{f, K) — > Jo x m\x\/Ji x ■ ■ ■ x M.[x\/ J^. 

Lemma 5.2. There is E P such that f = qo mod Jq. 

Proof. According to the argument presented above, V{Jo) Ci K = ^, hence 
there exists Uq E P such that —1=Uq mod Jq. This result is a special case 
of theorem 8.6 in |Lam] . 

We write / = /i — /2 for SOS polynomials /i = (/ + and /2 = 
+ i). Hence f = fi + Uo/2 mod Jq. Let go = /i + M0/2 e P- Then 
f = qo mod Jo- □ 

Lemma 5.3. f is a sum of squares modulo Ji, for all i = 1, . . . ,r — 1. 

Proof. According to the argument presented above, on each Wi,l < i < 
r — 1, / = > 0, and hence the polynomial u = f /ai — 1 vanishes on Wi. 
Then by theorem 12.21 there exists some integer k > 1 such that G J;. 
From the binomial identity, it follows that 

2 

i + u=\ > -(Tl^M 




The reader can see clearly in lemma 7.24 in |Lau] . 

Thus / = aj(u + 1) is a sum of squares modulo Jj. □ 
Now we continue the proof of theorem 15.11 



If Or- > 0, then by the proof of lemma 15. 3[ we imply that / is a sum of 
squares modulo J^. 

Lemma 5.4. If ar = 0, then there is qr E P such that f ^ qr mod Jr. 

Proof. By the assumption that / satisfies the BHC at each zero of f on K 
and by theorem 14.11 there exist g E P and h G M.[x] such that f = g + hf'^, 
i.e., f{l — hf) = g. Since / vanishes on Wr, f^ G Jp for some positive integer 

m—l 

m. Lett = hf,v = ^^t\ Thent,v G M[x],t"' G Jr, and (l-t)f = 1 mod Jr. 



i=0 

By the binomial theorem, there exist CjGQ, z = 0,l,...,m — 1, such that 

^m-l X 2 

^ Cif 1 mod Jr- 

■ i=0 ^ 
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This yields Qr & P satisfying 

/ = /(!- hf)v = gv = qr mod J^. 

□ 

To finish the proof of theorem I5.H we claim the following lemma. 

Lemma 5.5. Given qQ,qi, . . . ,qr G M.[x], there exists q G such that 

q — qi G Jj, Vi = 0, 1, . . . , r . Moreover, if each qi G P, then q G P. 

Proof. The proof is by induction on r > 1. Assume r = 1. As Jq + Ji = 
M[x], 1 = Uq + Ui for some Uq G Jo,Ui G Ji. Set q := Mq^i + ufqa; thus 
q E P. Moreover, q - q^ = ulqi + qo{ul - 1) = uggi - Uo(ui + l)go e Jo- 
Analogously, q — qi G Ji- Let t be the constructed polynomial, satisfying 
t — qo E Jo and t — qi E Ji. Consider now the ideals Jq H Ji, J2, . . . , Jr- As 
(Jo n Ji) + Jj = M[x](z > 2), we can apply the induction assumption and 
deduce the existence of g G M\x\ for which q — t E JonJi,q — qi E Ji{i > 2). 
Moreover, q E P if t, q2, qr E P, which concludes the proof. □ 

Using lemma 15. 2[ lemma 15. 3[ lemma 15.41 and lemma 15. 5[ we imply that 
there is g G P such that f = q mod /(/, K), i.e., f E P + I{f, K). □ 

Remark 5.1. If we replace the BHC condition by the radical condition of 
/(/, i^'), then we will have the following result. 

Theorem 5.2. Suppose that 

(i) f >0 on K, 
(a) I{f,K) is radical. 
Then f E P + I{f, K). 

Proof. From the proof of theorem 15. ![ by our definition of irreducibility, 
each Wi is conjugate symmetric (i.e., a point X G C" belong to Wi if and 
only if its complex conjugate X E Wi). By lemma 1 in |NDSj, there exist 
polynomials po.pi, . . . ,pr G M[x] such that PiiWj) = Sij, where 6ij is the 
Kronecker delta function. 

We consider the polynomial 

r 

i=l 

where go is as in lemma 15.21 By construction, q E P. 

Moreover, f — q vanishes on C{f,K), since f{x) = qo{x) = q{x) for 
X eWq (by lemma lOll and f{x) = ai = q{x) for X E Wi, Vi = 1, . . . , r. 
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By the assumption that /(/, K) is radical and using Hilbert's Nullstel- 
lensatz (see in |CLO] ). we deduce that f — q E I{f,K). This imphes that 



Remark 5.2. If we leave both the radical condition of /(/, K) and the BHC 
condition, then we will have the corresponding representations of strictly 
positive polynomials. 

Theorem 5.3. If f > on K, then f eP + I{f, K) . 

Proof. This follows similar argument in the proof of theorem 15.11 However, 
we can assume ai > ■ ■ ■ > a,. > 0. Thus, by lemma [573| / is a sum of squares 
modulo Jj, for alH = 1, . . . , r. Also by lemma and lemma [575| we imply 
that there is g G P such that f = q mod /(/, K), i.e., f e P + I{f, K). □ 



In this section, we present a result that is similar to theorem 4.1 in |DNP] 
and theorem 6.1 in |Hij . 

We consider the following optimization problem: Find 



In the case where K is compact, the SOS methods are based on represen- 
tations of positive polynomials on compact semi-algebraic sets, which were 
presented in the theorems of Schmiidgen |Schm] and Putinar |Puj . However, 
these theorems do not hold in the case where K is not compact. A more 
traditional approach in numerical optimization methods uses the first order 
optimality conditions. Using theorem 15. H and theorem 15. 3 [ we combine these 
two methods to give a procedure for approximating /* in the case where 
the semi-algebraic set is not necessarily compact. 

In order to implement membership in P + /(/, K) as a SDP, we need a 
bound on the degrees of the sums of squares involved. Thus, for (i G N, we 
define the truncated preordering as follows: 



/GP + /(/,K). 



□ 



6. Applications in optimization 



(1) 



r := inf /(x). 



Pd 




and the truncated critical ideal as follows: 
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Then we define a sequence {/^} of SOS relaxations of the optimization 
problem ([1]) as follows: 

(2) /rf* = maxr, 

(3) s.t.f{x)-TeP, + W,K). 

Obviously each T feasible in ([3]) is a lower bound of /*. So < f*. When 
we increase d, the feasible region defined by ([3]) is increasing, and hence the 
sequence of lower bounds {/J} is also monotonically increasing. Thus we 
have 

/r < /2 < /3 < • ■ ■ < /*• 

It can be shown that the sequence of lower bounds {/J} obtained from ([2]) 
and ([3]) converges to /* in ([1]), provided that /* is attained at one point 
X* G K. We summarize in the following theorem: 

Theorem 6.1. Assume f has a minimum f* := f{x*) at one point x* G K. 
Then lim^^-oo fd=f*- Furthermore, if f satisfies the BHC at each zero of 
f — f* in K , then there exists some G N such that f2 = f*, ^-e-, the SOS 
relaxations ^ and ^ converge in a finite number of steps. 

Proof. The proof is similar to that of theorem 6.1 in [Hi] (see also theorem 
4.1 in |DNP] ). However, we only consider the polynomial ring M\x\. □ 

Acknowledgment. The author would like to thank Prof. Murray Mar- 
shall, Prof. Ha Huy Vui and Assoc. Prof. Pham Tien Son for many inter- 
esting and helpful discussions on the topic of this work. 
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